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THEORY OF PLANE, SYMMETRICAL INLET PIFFUSERS* 
By Valter Bri'del 

PART I 


I - TFR TBREE -PARAMEDICAL GROUP OF THE INTAKE 

1 - Introduction. 


The present report ties in with the investigations on inlet diffusers 
by P. Ruden. The theory developed by Ruden (reference l) had produced 
results which found excellent confirmation in wind-tunnel tests (refer- 
ence 2) and in spite of certain still- existing defects, are technically 
very promising. The reasons for the new theory of the diffuser forms 
indicated by Ruden are twofold: first,' the arguments adduced in refer- 

ence 1 deal only with one specific operating condition, that is, a certain 
ratio of mean velocity within the diffuser to flying speed, while in the 
present report any desired velocity ratios are involved; second, a 
different choice of parameters and the increased possibilities of 
variations result in diffuser forms which cannot be reconciled at once 
with Ruden’ s theory. The first enables a theoretical check of the 
measurements made with Ruden 's diffusers at variable velocity ratio, the 
second permits the calculation of diffuser types which in many respects 
are superior to Ruden 's diffusers. 

So, while the present report seems to be a supplement and continuation 
of Ruden 's report (reference l), it is nevertheless a study by Itself 
and does not rest on the previous knowledge of Ruden’ s theory, with 
exception of the first part of section 6. The aerodynamic problem involved 
is the following (reference 1, page 3) : 

A certain volume of air from the air stream is to be intercepted 
and by conversion of speed to pressure or pressure into speed to be 
conveyed to certain airplane accessory devices. The manner in which the 

Wi 

ratio of the velocity inside the diffuser to flying speed — Is regulated 

Woo 

is not Involved here; the most important thing is to find the best possible 

*"Zur Theorie ebner, symmetrischer Fangdiffusoren. " Zentrale fiir 
wissenschaftlich.es Berichtswesen der Luftfahrtf orschung des Generalluf t- 
zeugmeisters (ZWB), Berlin-Adlershof , Forschungsbericht Nr. IV75/1 
and 2, September 20, 19^1* 
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shape for the contour of the nose. This, on the other hand, postulates 
that the usually inevitable increases of speed be as low as possible, and 
also that the wall thickness of the device be a minimum. It is found that 
these two requirements are in a certain contrast to each other; generally 
speaking, a reduction of the increase of speed at the diffusers in question 
can be attained only by an increase in wall thickness. One of the princi- 
pal tasks is to analyze this relationship. 

In conformity with reference 1, the following Idealizations 
are effected: The flow is assumed to be incompressible and perfect, 

so that it can be represented in the conventional manner by analytical 
functions of a complex variable; the diffuser is assumed to extend to 
infinity in flow direction. 


2 - Hydrodynamic Mapping And Introduction of the 
Parallel Strip as Reference Area 


Figure 1 represents the contours of a plane, symmetrical inlet 
diffuser and at the same time indicates the qualitative performance of 
a corresponding flow. The problem is to change over from' this schematic 
figure to a quantitative one. This is accomplished by the method of 
hydrodynamic mapping. 

Hydrodynamic mapping is, as Is known, conformal mapping in which 
one space, along with the flow In that space. Is transformed into 
another. Three functions are involved in the analytical representation: 
The stream function (the complex stream potential) F(z) of the original 
space, the mapping function £(z) which Is considered to take a space 
on the z -plane to a space on the £ -plane , and the stream function $(£) 
of the transformed region; F(z) =$(£)* that is, the values of the 
stream function are simply transplanted from the z- to the t, -plane . 

The "complex velocities" w(z) and cd(£) are given by w = ~ and 

dz 

cd = ~ and related accordingly by 
dQ 


w dz - cd d£ 


or 


( 1 ) 



w 
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It is to be borne in mind that only the conjugate complex value of 
the "complex velocity" represents the velocity vector. 

A flow defined by an analytical function F(z) can be mapped 
hydrodynamically in a number of ways. The two most frequently employed 
mapping methods are : 

(1) £ = F(z), in which case $ (£ ) = £ , that is, the mapped space 
simply contains the basic flow, a horizontally directed parallel flow 
with the constant velocity 1. Any flow can thus be transformed to the 
basic flow. 

(2) 5 = Msi = w(z). This is the holograph mapping method 

dz 

successfully applied by Euden to the inlet diffuser problem. 

In the following, a different mapping method is used, since, as 
pointed out in the introduction, an infinite series of stream functions 
corresponding to the different velocity ratios wi/wc must be considered 

for each' inlet diffuser, rather than a single stream function, with 
which, of course, the above special mappings are invalidated. The image 
space is so chosen that the group of stream functions permits the 
simplest .possible representation. A few trials confirm that the parallel 
strip affords the simplest solution. 

Assume that the function £ (z) maps the simply connected space 
situated outside the contour of the inlet diffuser single-valued and 
conformally on to a parallel strip bounded by the straight line T(£ ) = ±it 
The line of symmetry of the z -space, taken parallel to the real axis, 
is to become the real axis of the £ -plane, and the direction 
of flow is to remain the same. In the z -plane the line 
of symmetry connects the two infinitely remote boundary points of the 
space; they become with their surroundings the two tips of the parallel 
strip. Figure 2 represents the flow in the £ -plane diagrammatically; 
the general function $ (£ ) which produces such a flow is. given by 

®(£) = Ql£ - c 2 e ”^ (c-l > 0, c 2 > 0) (2) 


To find this expression, the -parallel strip is mapped by =• 

on the plane cut along the negative half of the real axis 
and the flows arising in the £ -plane calculated. Since $ '(£) by 
analytical continuation beyond the strip edges obviously has the 
period 2jti, the derivation of the transplanted stream function 

must be unique. Its singularities must be looked for at 0 and w; 
at 0, corresponding to the left strip .tip, a. source and dipole appear 
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combined, while at 00 arising from the right tip, a sink occurs. As a 
result, the general expression of <I>i( £j_) can he immediately expressed 

in the form 3>i(£) = a log £l + -p-, with positive a. The outflow 

*1 

direction of the dipole points to the right, that is, b must be negative 
(real) . Putting a = ci, b = -c 2 and reverting to the £ -plane, gives 
exactly the above expression. 

The stagnation points are obtained by setting the derivative 
equal to zero. Prom 


*'(£) 


ci + c 2 e 


-£ 


= 0 


follows e~ £ = Putting £ = I + irj and accordingly 

>“£ = e”^ - ^-7 = 0 “^ 


e = = e 
e 3 cos tj •= 


-6 ci _ | 

=» r>na r\ - i. ^ = 

C 2 


e“ 5 (cos rj - i sin tj) gives the conditions 

t 

and e“ sin tj = 0 for the stagnation points. 


Therefore, sin tj = 0, while cos tj is negative, which in 
limiting t, to the parallel strip means tj = +n. Accordingly 

also: e" ^ = fi and I = log — (natural logarithm). On putting 

c 2 C 1 

£ = log + iti, that is, e^° = the expression for $ (£ ) can be 


written in the form 


C 1 


®(5) 




( 3 ) 


or since no additive constant is involved, 


«(5) = c i (5 -5 0 ) 



W 


The system of the streamlines is not ehanged by a change in the 
constant Ci; if £ 0 is varied - obviously T(fj Q ) must remain 

equal to it - the system of streamlines experiences a simple parallel 
displacement, as shown by the last equation. 
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3 - Determination of 1}he Differential Quotient of the 
Function z(£) Which Maps the Parallel Strip 
on the Outside of the Inlet Diffuser 


Given the function z(£) which maps the parallel strip on the 
outside of the Inlet diffuser, the transplanting of the general stream 
function $(£) makes the general stream function F(z) available. As 
regards the derivation of the function z(£) the following may he stated: 


Along the real amis, ^ is, on the whole, real and 
positive. Running through the straight line T(£) = it from left to 

right, assumes at first negative (real) values in correspondence 

with the straight part of the contour, then. In correspondence with 
the curved nose .contour, it assumes values which geometrically expressed 
belong to the lower half -plane, and lastly the values become positive 

(real) : Letting t shift leftward to infinity, must, in amount, 

d£ 

increase beyond all limits, while at unlimited distance toward the 


right, 


tends toward a finite (positive) limiting value. 


Further 


identification is not possible without a certain arbitrariness. The 


7 ; y 

following assumption is made. Suppose quantity when t, passes 

through the upper right edge of the parallel Btrip, runB through a 
line whose center piece located in the lower halfplane is half of an 
ellipse symmetrical to the real axis (fig. 3 ) • 

Altogether, considered as function in the parallel strip, gives 

then a conformal mapping on the shaded space represented in figure 3- 
For the first, the mapping of the total boundary follows frcm reasons 
of symmetry, then, the Inside mapping according to the principle- 
of the characteristics of the boundary. 


Of the multitude of potential diffuser forms, one series defined 
by a f ini te number of parameters -was selected. The ellipse is, for 
shape and position, given by three real parameters^ the right-handed 


end point of the line, that is, 11m gives a fourth parameter 

and a fifth is ultimately afforded by the possibility of effecting 
any desired parallel displacement in the £ -plane In direction of the 
read axis. Thus five real constants enter the general 


dz 

dr 


But, since the parallel displacements in the £ -plane 


function 
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are unimportant and M Itself may be multiplied by any positive factor 

without modifying the respective diffuser form, the factual result 
is a three-parametric continuum of inlet diffusers. 

The analytical expression of is obtained by mapping the 

d£ 

parallel strip through 

t = e“£ (5) 


on the t -plane cut along the negative half of the real. axis. 

The right tip takes the vicinity of the zero point, the left tip that 
of the infinitely remote point. The next problem of transplanting the 
plane thus cut on the space of figure 3 is essentially synonymous 
with the problem of mapping the plane fitted with a finite straight 
slit on the outside of an ellipse. The solution is predicated on the 

knowledge of the napping attained by z + — . This function transforms 

2 

the outside of the unit circle in the plane cut along a finite 
strip, and it maps the outside of the circle | z | = R > 1 on the outside 
of the ellipse; the infinitely distant point remains fixed. The outside 
of the circle | z | = r < 1 changes into an elliptically bounded space; 
this, however, lies in part as two-lobed space above the plane, which 
is to be avoided in the present instance. To pass from the slit area 
to the outside of an ellipse involves essentially, that is, apart from 
similarity mapping, the application of the inverse function of z + i 

which gives a circular space, then, after a second similarity mapping, 
the transfer to the elliptical boundary according to z + i. Through 

the similarity mapping certain real constants enter in the formulas, 
which must be subjected to certain restrictions in order to be certain 
to obtain an, on the whole, single-lobed image space. The mode of 
calculation is as follows: 

Assume that the end points lie at -t^ 
that case the slit is widened out 'by 


t' = t + y(t + t]_)(t + ± 2 ) 


and -tg. (t 2 > t-^ > 0). In 


+ 1 2 


2 
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in a circle the center of which is the zero point of the t' -plane. This 
result is easily checked, "because for -tr, < t < -t^_ the square 
root is purely imaginary, and 


1 1 ' | = 1 1 + ~ - + ^(t + t-^ (t + t 2 ) 


" (t + t-j_) (t + t 2 ) 


[t + 


tl + tr 


" t l t 2 “ 


t 2 - t-|_ 


that is, jt'J = constant. The function t" = at ' + with a > 0 

and h = 0, obtained from z + — by similarity mapping results in the 

general ellipse whose axes coincide with the coordinate axes. To the 
right-hand end point of the diameter of the circle, hence 

to t* = — — , there must correspond a "positive" value, that is. 


a”— — + — > 0 (6) 

2 ta - ti 

in order to prevent a two-lobed overlap of the plane at any point by tbe 
image space. Likewise, the uppermost point of the circle, 

t' = 1 — must take a point on the upper half plane, that is, 



2b 

t2 - ti 


> 0 


a 


2 


(7) 
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The prefix of b decides which of the 
principal axis. The desired quantity 

to t " + c (c real) . 


coordinate axes becomes the 

most generally be put equal 



t + 


t-L + t 2 


+ 


\j(t + tq)(t + 



t + ; , i *2 + \J( t + tq)(t + ,to) 


+ c 


= a t + 


P \j ( t + t-^Ct + t2) + 7 


(8) 

> 


with 


a = a + 


4b 

(t2 - tq)? 


tl + to . 

7 = a — - - + c > 

2 • 


P = a - 


4 b 


(t 2 " t x ) 2 


( 9 ) 


The expression of 


dz 

dT 


therefore actually contains five real 


constants, which must satisfy the following conditions; first, 
tg > ti > 0 . The conditions (6) and ( 7 ) are equivalent to a > 


P > 0; a = °? * P itself then becomes positive. Lastly, the 

ellipse must comprise the zero point even after the tr ansla tion 
effected by the constant c. This is the case only when the values 

dz 

of ^ relating to t = -tj_ and t = -tg are positive and negative, 
respectively. From 


-a t-j_ + 7 > 0 
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and 


■a + y < 0 


follows 


a < 7 < a tg 


( 10 ) 


The square root must be calculated positive for positive t. In the 
analytical continuation along the dissected t -plane, the values 
obtained for t < -tg are negative, while the values in the 
interval -t2 < t < -t]_ relate to the positive and negative half plane, 
depending upon whether the real axis is approached from above or below. 


4 .- Determination of the Mapping Function Itself by 
Integration. Calculation of the 
Contours of the Inlet Diffuser 

The function z(£) Itself is found by simple integration. From 

= at + P \j ( t + t x )(t + tg) + 7 

and t = e =, that Is, t = - log t follows 

dz _ dz = 1 dz 

dt d£ dt t d£ 

and 


-a 


- 1 


t t 


■+ \j( t + t]_)(t + ±*2) 


\J { t 4 - tl)(t + t 2 ) ^ 


z = -at - 7 log t 
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The last Integral is changed, "by the introduction of 


f t + ijj 

t + tr 


, as auxiliary variable to an integral of a rational function and 
computed direct. Altogether the result is 


z = - at 


- 0 \j( t + t x ) ( t + tg) - ( 7 + 0 \/tlt 2 j log t 


0(ti + t 2 ) log ^ \Jt + ti + + 20 \Jt ± ± 2 log (yt 2 (t + tj_) 

(li: 

* 

^t^(t + tg) j + constant 


The square roots containing t are all chosen positive for 
positive t -values, with which they are then unequivocally defined over 
the dissected t-plane in the sense of the analytical continuation; 

\J tqt 2 is positive also. For the logarithmic functions which, like the 

square roots in the dissected t-plane, are unbranched, the principal 
values, that is, the values whose imaginary parts lie between -jti 
and jti, are visualized as being inserted. If the additive constant 
is chosen real, the axis of the inlet diffuser coincides with the 
real axis. 


On letting t run along the lower rim of the negative 
half of the real axis, z describes the upper contour of the desired 
diffuser. The Intervals (0, -tq) and (-tg, -«} give the straight 

sides, the center piece (-ti, -tg), the nose of the diffuser. 

Putting z = x + iy. 


x = 


-at - 7 log 1 1 | -+ constant, 

0 \/(-t-]_ - t)(t 2 + t) + 0 (tq + tg) arc tan 


y-fa 1 - 1 

\[t2 + t 


20 \J t-j_tg arc tan 


VtgCtj^ - 1) 

^(tg + t) 


+ constant 
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Is valid along the center piece , or, if t = -r (t-]_ = t = ± 2 ) and "both 
additive constants are taken at zero 


x = aT - 7 log T 


y = 


P \fl t - ti) (t2 - t) + + ±2) arc tan ^ 


\J*2 " T 


-P\/tlt2 arc ^ arL 


JSHEI 

X/tiCtg - T) 


( 12 ) 


The roots must "be extracted positive, and if the argument, as here, 
varies between . 0 and °°, the principal branch from 0 to ^ is 

chosen for the arc tan functions . 

The formulas (12) define the upper half of the diffuser. 
Horizontal straight lines are drawn toward the right to infinity from 
the two end points of the nose given by t = t j_ ' and T = tg. 

x = at-L - 7 log t^_, y = 0 


and 


x = atg - 7 log t 2 , y = P(t x + t 2 - 2 \Jt 3 t 2 ) | 

The complete determination of the diffuser must include the inside 
width 2h; from (ll) follows 


h = (7 + P tit 2 ) it 


Of the five constants which define the individual diffuser only 

7 

three have any essential significance. If ==■ is multiplied by a 

d b 

constant positive factor and a real constant is added to £ , the 
diffuser form remains the same. But adding a real constant to £ is 



12 


KA.CA TM 1267 


reflected in t as multiplication "by a positive constant. So, without 
causing any substantial change, the expression 


at + 3 \J ( t + ti)(t + tg) + 7 (13) 

can be multiplied by X > 0 and t itself be replaced by |at(n > 0) . 
The expression (13) then becomes 

aX|at + pXp, \j(t + ii^t^Kt + n^^tg) + ^7 (l4) 

Since a and 7 are positive, because of (6) and (10), X = — 

and n = ■ 21 and (l4) reduces to 
a 

t + 1 + 3 ' \yf^7^7(tTvy (15) 


with 


3' = K t* = S t lf 
a x 7 



or expressed in different notation 


= t + 1 + 3\/(t + t x )(t + t 2 ) 

The three quantities 3, t-j_, and t2 satisfy the conditions 

tg > 1 > ± 2 . > 0 


and 


(16) 


3 > 0 
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So the formulas for computing the contours read 
x = t - log T 


y = 0 \J( t - t-j_) (tg - t) + 0(t-j_ + tg) arc tan — - — I— ^ 

n/^ ' T 

- 2 P\/t t arc tan " — (17) 

/12 Vtl(t 2 - T) 

h = (1 + Px/tjt^) it 

5 - Calculation and Representation of the Telocity 
Distribution Along the Contours 

The general inlet diffuser flow is obtained by transplanting the 
general function 0 (£ ) = C]£ - c^e - ^ according to (5) and (ll) in the 

z -plane. The magnitude of the velocity along the contours is the 
principal point of interest. The complex velocity in the z -plane is 
given by 


v(z) = §* = 

dz dT d£ 


(18) 


where 


= °1 + °2 e ~^ = c i + c 2^ arii ^ % = t + 1 + 3 \j( t + t-^) (t + tg") 


d£ J- J- * — dT 

in the simplified method of writing. At the boundary t is negative 


(real) . Draw a t-axis and plot the amounts of 


as ord in ates. The graph of 


d® 


d® 

d.5 


and against it 


consists of two straight lines, 


or more exactly, pieces of straight lines, rising from the 

point t = — = of the t-axis at the same angle toward both sides; 
°2 
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one straight line -withdraws toward the upper left to infinity, 
terminates at height C]_ of the ordinate axis. The amount of 

can he formularized as follows: 


the other 
dz 


d£ 


t + 1 + 3 \/( t + tx)(t + tq) for t > ti 


-< 


\/(t + l)^ - p 2 (t + ti)(t + tg) for -t]_ ^ t sT -t 2 


( 19 ) 


-t - 1 + ( 3 \/(t + tj_)(t + t 2 ) for t g -t 2 


All square roots are positive. Accordingly 


dz 


is represented by 


a curve that appears built up from three arcs. The two outside ones 
are hyperbolic arcs, while the inside one can relate to any conical 
section and can even be rectilinear. 


Individually the following applies: Conformably to the 
condition t 2 > 1 > tj_ > 0, tq_ and t 2 are chosen fixed, so that 

the two end points of the central arc which are at the same time the 
connecting points of the outer arcs are defined. The ordinates 
are 1 - t^, and t2 - 1 . On examination of the entire group of 

curves obtained for variable 3 (fig. 5) it is seen that the case 3 = 0 
regarded as limiting case, since 3 must be positive by assumption, 
results in a pair of straight lines. As 3 increases the ordinates 
along the entire line increase. The outer hyperbolic arcs deviate in 
their connecting points with vertical tangent and become steeper and 
steeper with increasing 3 . The inside piece cancels out hyperbollcally 
which is, the' t-axis is the secondary axis for the first appearing 
hyperbolas. The curvature of the arc decreased continuously . If 3 
reaches the value 


Po = 


gy/U - ti)(tg - 1) 


t 2 - tj_ 


(20) 
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the hyperbolic arc becomes the straight connection of the end points; 
at further increasing (3 the arc becomes concave downward. In 
general, 0 O is less than unity. The 0 values lying between 0 O 

and unity give hyperbolas for which the t-axis is principal axis . 

0=1 gives a parabolic arc with the t-axis as axis of symmetry, 
and 0 > 1 results in elliptic arcs with progressively increasing 
steepness. The special case 0 O = 1 occurs only for + tg =2. 

The end points of the middle arc are then at the same level and the 
convex hyperbolic arcs appear ing for small 3 0 pass over the straight 

line corresponding to 0=1 directly in .concave elliptic arcs . 

The amount | w j of the boundary velocity follows by (l8) by 
division of the two functions of t; the numerator function is 
represented by the straight line, the denominator function is represented 
by the triple product of conic -section arcs. 


6 - Comparison With Ruden' s Investigations 

Ruden (reference 1) used the holographic method of mapping throughout 
his experiments; the most general space of the w -plane (w = complex 
velocity) taken into consideration by him is a circle the center of 
which lies on the real axis and which exhibits two 
radial incisions along the axis. The space is shown in figure 6. 

The slot ends w*, and w-j_ correspond to the flight speed and the 
terminal speed inside the diffuser, that is, the speeds relating 
to t = 00 and t = 0. The values w-^ and w 2 indicate the velocity 

at the end points of the curved part of the diffuser wall. According 
to Ruden: 


w 2 = > w ± > 0 = wi (21) 

The analytical relation between w and t is readily indicated. 

The discussion is restricted to the lower half of the t-plane and of 
the w-circle, which are clearly and conformally referred to each other. 
Excluding, in the first instance, the appearance of equality signs 
in (21) gives the point coordination 


t 

0 “ 

_t l 

-t 2 

w 

Vi w„ 

wi 

Wp 


tj_ and ±2 > ti signify any two positive values. 
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Mapping the lover half plane by t 1 = — . + ^ in itself in such 

t + t2 

a way that -tj_ and -tg change to 0 and <», vhile the 7alue 

t' = 1 corresponds to the value t = 0 the formation of the square 

root t " -A 7 results in a quadrant (main branch of the square root 
for positive t') and a linear mapping with real coefficients, vhioh 
are written in the form, 


w 


= at" - b 
ct" + d 


finally gives the desired w-space 


Accordingly 


a \Jt + t-j_ - bsjt + 

c \J t + t]_ + d\/ t + tp 


The above table of values then appears as 


( 22 ) 



0 

\/ti - \[^2 

\J iq_ + d \J^2 


00 -tj_ 


a - b _ b 
c + d d 


-t 2 


a 

c 


The value w = a> corresponds to a negative t". From it and from (21) 
the following conditions for the coefficients a, b, c, d are deduced: 


-> 0, y>0, -> 0 

c ’ d * c 


a \[t[ - b \[t 2 

’ c^I + dV^ 


> 0 


( 23 ) 
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The first three inequalities indicate that all coefficients have the 
same signj without restriction of generality they shall "be positive. 
The last inequality signifies then 



{2k) 


By making the numerator in (22) rational, 


(a 2 - b 2 ) t + (a 2 tq - b 2 tg) 

w = — — - 

(ac + bd) t + (ad + bc)\/(t + t]_)(t + t^ + (acti + bdt 2 ) 


Quantity a 2 >- b 2 is positive by reason of (2l+) . Since the 
coefficients a, b, c, and d may be multiplied by any positive 
constant, without modifying w(t), 8.2 - b 2 = 1 (25) he assumed. 

Therefore : 


with ■ 


t + t. 


w = 


At 


+ B \f { t + t]_) (t + t 2 ) + C 


(26) 


t Q = ~ h 2 t 2 

A = ac + bd 

(27) 

B = ad + he 

C = act]_ + bdt 2 



18 


MCA TM 1267 


A, B, C, and t c are positive. On comparing the expression for w in 
(26) with 


Cgt + Cj_ 

at + 0\J(t + t 1 )(t + tg) + y 


which can also he written in. the form 


t + t r 


w = 


c 2 


at + P\/(t + t 1 ) ( t + tg) + y 


(28) 


where the denominator is used again in full generality, it is seen that 
all the forms considered by Ruden are included. The only conditions 
which the present coefficients, aside from that of being positive, must 
satisfy are: ^2 > ^1 > 0 and (10); tg > t^ > 0 is also satisfied 

for (26) and from (27) follows immediately the inequality corresponding 
to (10) 


Atj, ^ C *C At g 


Two more facts stand out: 


(1) The numerator in (26) is defined when the denominator is 
known; this is proved by (25) and (27). But the zero place of the 
denominator in (28) remains arbitrary for given numerator (<0). 

(The factor c 2 is naturally unessential.) The limitation in (26) 

is associated with the nature of the hodographic method, which momentarily 
comprises omy a specific operating condition and does not detract from 
the generality as far as the diffuser forms are concerned. 

(2) It is readily seen that for specified values of tp and tg 
the formulas (25) and (27) do not permit the ratios* A : B : C to vary 
in the same manner as the ratios a : 3 : 7. Thus at tp = 1 and tg = 2 

v 

the proportional equation A : B : C = 1 : 1 : “ can certainly not 

be -made to agree with (27), while being able to put a = (B = 1 and 

y - -2.. To this extent the formula is also more general as regards 
2 
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diffuser forms than Ruden's formula (in reference l) although, the 
number of constants Is not greater; Ruden's study of diffusers showing 
a constant speed along the nose contour In normal operating condition 
requires, according to the present theory, linearity of the curve 

of j representing the central piece . The normal operating conditions 

further require that the zero place of the numerator function of (28) 
lie, geometrically speaking, in the straight extension of the central 
piece. The diffusers characterized in the foregoing by 0 = 0 O correspond 
accordingly to these special inlet diffusers. The degree of generality 
of these investigations surpassing Ruden (reference l) in two ways, is 
especially clear. On the one hand, the numerator function in the 
expression of w need not disappear exactly at the t-place lying in 
the extension of the center piece (more general operating condition), on 
the other, this extension need not as in Ruden's report, intersect the 
t-axis at all in a point with negative abscissa (more general diffuser 
forms ) . 

II - THEORY OF OPTIMUM DIFFUSERS 
1 - The Characteristic Quantities w-j_, w^, and v*5 


Earlier in the present report the complex velocity of the general 
inlet diffuser flow was defined and the amount of the velocity along 
the diffuser contour analyzed. This amount is hereinafter designated 
by wj the complex velocity previously denoted by w does no longer 
appear. With s(t) as the magnitude of the quantity ~ along the 
boundary ' = 

t + 1 + 0\/(t + t;i_)(t + tg) for t > -t-L 


s (t) = <^\/(t + l) 2 - P 2 (t + t-L)(t + t2) for -t x > t ^ -t 2 (29) 

-t - 1 + 0\/( t + t 1 )(t + t 2 ) for t < -t 2 


It is borne in mind that t assumes only negative values, 
as function of t 


So far w 


. . Ci + Cot 

v(t) = 1 , 

s(t) 


t + t r 
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In the subsequent study only velocity ratios will "be the controlling 
factors. Without restriction of generality C 2 can he put as 02 = 1. 

To the values t = 0 and t = -» correspond the velocities and 
they are readily obtained when t Q , tq, tp, and (B are known. On the 
other hand, it is quite difficult to set up a general expression for 
maximum increase of speed w^-y relating to the given parameter values. 

In each specific case wmax is also easily calculated. To establish the 
abscissa value for which w attains its maximum, consider a pair of 
straight lines (fig. 8), radiating from point (-t Q , 0) and symmetrical to 
the vertical t = -t Q . Assume that this set of straight lines, starting 
from a very flat slope with respect to the t-axis, becomes progressively 
steeper until it meets the curve C representing s(t) for the first 
time. The intersection point gives the desired abscissa value. Actually, 
the pair of straight lines can be regarded as geometric^ representation 
of the numerator of w in the particular limiting position, since no 
proportionality factor is involved. In that case w < 1 and the equal 
sign is. reached exactly in that particular point. Naturally, can 

occur at several points. The difficulty of a general determination 
of VmfiT arises from the faot that different cases have to be distinguished, 
depending upon whether the maximum point lies within the central arc 
of C or at a point of discontinuity. The general examination is 
restricted to the determination in place of v TOT to the higher of the 
two velocities which correspond to the two points of discontinuity, 
indicated as w* and noting that in many cases w* = w ^^ . For the 
equation w* = w mRy certainly applies when C has a straight or even 
slightly convex central piece, and indeed, it is then fulfilled for 
all values of t Q . It also applies for slightly concave central piece 
in certain conditions. The subsequently discussed optimum diffuser 
forms show this characteristic practically generally, which physically 
implies that the maximum increase of speed occurs at one of the two 
end points of the diffuser nose. At any rate the difference of w mflx 
and w* in the diffuser types treated in section IV is significant. 

The value of t = -t^ or t = -tg corresponding to the higher speed is 
readily - apparent . At fixed s(t) and variable t$ a reversal takes 
place when t 0 passes through unity, or when the apex. of the two straight 
lines falls In the straight extension of the line connecting the points 
of discontinuity of s(t). With -t^ signifying the abscissa of the 

point of Intersection of this extension with the t-axis 

t = St^tg - tq - t g 
3 tj_ + t2 - 2 


( 31 ) 
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After ex amina tion of the different positions successively the 
result is as follows: 


(a) t 3 < 0 

(a) t Q % 1 v ± : w^ 

(P) t Q > ! v ± : v a 

(h) 0 ^ t 3 < 1 

(a) t Q -g t 3 v ± 

( P) t 3 -g to < 1 w ± 

(7) t Q > 1 

(c) t 3 > 1 

(a) t 0 ll w ± 

(P) 1< t Q t 3 w ± 

(7) t Q > t 3 w ± 


V* = 


1 + 


p\/t^ 


. 1 . t g " % 

■ 1 + P ■ t2 - 1 


wv = 


t o . 1 . t o ~ t l 


1 + p^t„t. 


1 + p 1 - tn 


12 


( 32 ) 


w « : w * = 


Woo = = 


Woo : w* = 


w : v* = 

00 * 


v : w* = 

00 


V co *• W* = 


1 + p^ 
*o 

^ 1 + 13 1 - t l. 

1 . * 2 .; s o 

1 + 

*o 

.. . 1 . V- l i 

1 + p \l 
% 

^ ‘ 1 + P ' 1 ' *1 
1 . *2 “ *0 

1 + 

% 

, r- * 1 + P * t p - 1 

tlt 2 

, t_ - t, 

l.o 1 

i + e\j 

*o 

^ 1 
1 . *0 " t 2 

1 + pV^ 1 1 13 ~ 1 


( 32 ) 
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The case t^ = 1 does not occur; t^ = °°, that is tq + t 2 = 2 
can he incorporated under (a) or (c) (a, (3) . 


The ratio of the speeds Wq, w w , and w* permits an appraisal 
of the mode of operation of the inlet diffuser. In addition the 
wall thickness d, measured, say in ratio of l/2 inside width, is 
important. By (ll) 



( 33 ) 


The question posed in the introduction concerning the best possible 
design of diffuser contour can now be formulated as follows : 

The parameters t^, t 2 , and 3 shall be so defined that 

(1) The quantity 5 becomes as small as possible 

(2) At variation of Wq : w^ within a specific range, that 

is, at variation of to within a corresponding interval, the quantity 
Wjuax (or w ) is the smallest possible with respect to the higher of 
the two velocities w M and vq. 


2 - The Question of Most Favorable Choice of Parameter 

' A mathematically precise treatment of the optimum problem presents 
great difficulties even with w m .r instead of w*. There are two 
reasons for these complications. First, it is, of course, easy, for a 
fixed ratio v. : w M , to make the problem a precise minimum problem. 
Either assume fixed 5 and define tq, t 2 , 3, and t Q in such a way 

that (33) is complied with and — becomes a minimum, or else specify 

w_, 

T J 00 

* 

“ and attempt to reduce 5 to a minimum. But at variable w-» : w„ 

oo v* 

on the other hand, — is naturally variable also, and a precise 

formulation of the minimum problem is not possible without option. 
Second, the complicated form of the function s(t) pieced together 
from several analytical functions is disturbingly noticeable in the 
calculations and the attempt to restore the organic character of a 
uniform function introduces new difficulties. For this reason it was 
decided to set up rules based upon the geometrical view according to 
which the parameter values are chosen to suit the purpose. Everything 
else is left to the special numerical calculation. 
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(a) The case of constant velocity ratio : w^ is analyzed first. 

The numerator function of the expression for w, represented geometrically 
hy a pair of straight lines (g 1 ' J g^'), is visualized as fixed. If 

a certain increase of velocity is admitted, the curve C representing 
the den omin ator function s(t) must "be situated in a part of the (t, s) 
plane which is downwardly "bounded by a certain pair of straight lines 
(Si, &2^ • This P a Tr iLas the same vertex point as g 2 *) and- 

like (g]_‘, gg'), a vertical line of symmetry. Even (g-^, gg) can be 

regarded as geometrical pattern of the numerator function, since constant 
factors are not involved. The best diffuser is obtained when the 
infinite angular space available for C is utilized to the best advantage 
that is, the curve is as close as possible to the boundary. There are 
two possibilities. The first consists in choosing the central piece 
of C straight and placing it immediately on g^ or gg. This way 

a two -parametric continuum is separated from the three -parametric diffuser 
continuum defined by t-^, tg, and 0 ; 0 O is obtained ( 20 ) through tp 

ana t 2 . This two-parametric continuum decomposes in two partial 
groups j one, characterized by 82 as carrier of the straight part of C, 
is the group analyzed by Euden particularly. The two parameters t-^ 

ana t 2 permit, as will be seen, the realization of any velocity ratio 
w ± : w M ; Wjug^. For Euden' s group, w A < w*,. 

The second possibility consists in placing a point of discontuity 
of C on both g x and gg and choosing 0 very small so that the central 
piece of C clings very closely to g-j_ and g 2 . The limiting case 0 = 0 
leads direct to the pair of straight lines itself) the corresponding 
diffuser consists of two infinitely thin straight walls and represents 
the optimum solution for w-j_ = 

(b) Suppose the ratio — is variable. The case of unlimited 

Woo 

variation serves as basis, but wj_ and w m themselves are visualized 

as varying only between zero and a finite limit. In this event the 
fixed numerator function previously represented by (g-j_, g 2 ) is replaced 

by a two-parametric continuum of numerator functions. What part 
of the plane do the corresponding pairs of straight lines cover? 

The Int ersection points with the ordinate axis cover a finite 
interval starting from the- zero point, while the vertex points cover 
the entire negative half of the t-axis. To vertex points of very 
great distance correspond very flat pairs of straight lines, that is, 
very httw.i i values of . The measure of rise varies altogether between 

zero and a finite value as exemplified in figure 9 . The space filled 
by the pairs of straight lines is itself bounded by a pair of straight 
lines which again is denoted by (gi, g 2 ). 
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3 - Theory of Best Diffusers for Fixed w^ : w^, 

The diffusers for which the curve C exhibits a straight central 
piece are characterized by the condition 



2\J(1 - t-jHtg - 1) 

t 2 - ti 


(34) 


The line connecting the points of discontinuity meets the abscissa axis 
in point -t 3 , where t 3 shall now be positive. Putting t Q = t^ 

(31)) gives the normal operating condition of the diffuser, for the 
point (-t-^, 0) corresponds to the vertex point of (gq, g 2 ). Since 

w max = ^ 


V * TaT •XJ' 

• w » • "max 


2tit g ~ tp - t g 


(t-L + tg - 2)(1 + P^Vg) 


w i :w oo :v max = 


ti + tg - 2tqt2 


(2 - ^ - tg) (l + $\J tqtg) 


1 

1 + 0 


t 2 ~ H 

ti + tg - 2 


for t 3 


<1 


1 

1 + 0 


t 2 ~ t l 

2 - tq - tg 


for 


(35) 

t 3 >1 


for 0 the value (34) is inserted. 

Equations (34) and (35) ar e now used to compute tq, tg, and 0 
for prescribed ratios w^ : w M : w^^. From 


0 2 = 


4(tf + tg - 1 - t x t 2 ) 
(t2 " 4q) 2 . 


follows 



h : 2)g , (1 + 3)2 JkL 

(tg - tl) 2 VnR.y 2 
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hence 


or 


w. 


CO 


W, 


mas 


1 - P 2 = 1 - P 
(1 + P) 2 1 + P 


v. 


P = 


mas 


w . 


w, 


mas 


+ 


The application of this value gives 


(ti + tg - 2) g = W^ 2 v^ 2 

/ , , \2 / 2 2\2 
(t 2 - t x ) (v^ + v* ) 


( 36 ) 


or 


*1 * *2 - 2 


^max v °° 

* p p 

wmas + w oo 


(37) 


Let t 3 he less than 1, that is, + t 2 - 2 he greater than 0. In 

this case (37) carries the plus sign; and. (37) is written In the form 


hence 


( tg - 1) - (1 - t X ) _ gffmp.Y Vcc 

(tg - 1) + (1 - ti) ' 


t_ - 1 (w + v) 2 
2 _ _ mas ” 

1 “ tj_ ( ffimr “ Wco)2 
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According to this relation put 


1 - t x = XCw^ - wj 2 ; t g - 1 = + v^) 2 


or 


ti = 1 - - O 2 ; t 2 = 1 + + wj‘ 


Quantity w^ must he used for determining X. 


2t 1 tg - t x - t 2 _ 1 _ X ( v max 2 ~ v qq 2 ) 2 
tj_ + tg “2 2 V WELX w °° 


and after a simple calculation 


V 1 


2w, 


max 


Woo " ^ (^] nwy 2 “ v oo ) 


2\2 


w. 


max 


w 


max 


+ w„ 


(v 2 - V 2 )l 
vw max 00 ' 


1 + X ( W-rpyj/y- + w M ) 


1 " " O 


This equation contains only X as unknown. The removal of denominator 
and square root leaves an ordinary quadratic equation for X. Its two 
roots are 


^max^max 2 “ ^f^oo) (vco - w^) 
( w max 2 ” w oo 2 ) 2 (vmax 2 + w ± 2 ) 
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and 


2w_ 


Xg — 


(■»; 


V 

CO 


max 


22 

- V ) 


The second solution proves useless or at least dispensable . It 
leads to t Q = 0, a result that applies only to = 0 . But 

for w.j_ = 0, X^ gives a value that exactly agrees with X 2 . 
Therefore put X = X^, so that 


w maoc 2 + t *° 2 (wmax + ^i) 2 

’1 " 2 2 2 (38) 

w max + v i ( v rrw.-r + w M ) 


and 


v max 2 + v oo g ( v max ~ v j) 2 

2 2/ \2 
'Wmax + ^ w max “ 


For the present the calculation gives t]_ and t 2 in the form. 


tl = 1 


2W TTW.X^ V TTW.'ir ~ V l V Q3^ V oo ~ V j) 

2 2 2 
( w max + w i )( w max + w«,) 


to — 1 + 


2w (w 

mo nr * y 


2 _ 


max max 


w.w )(w 

1 00 a 


w i ) 


(39) 


(■w, 


max 


+ Wj_ ) ( w^y - Woo) 


These expressions indicate that t 2 > 1 > t-j_ only when w^ < w^; 

therefore (38) applies only to this case. This limitation is linked 
with the previously made limiting assumption t 3 < 1. As is imme diately 
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apparent from (39), + t 2 - 2 is actually > 0, that is, tg < 1 

for wi < w^,. For t-j_ + t 2 - 2<0 the minus sign is carried in (37). 

Then 

t2 - 1 _ (yarn: ~ W») 2 

1 ~ tp (wmajc + Vod )2 

and a similar equation as above gives 

Wmax 2 + v °° 2 t y max ~ v i )g 

"^1 — 2 2 / \ 2 

vmajc + *i (v MT - w*) 

(4o) 

v max + v °° ( v iaax + v j) 

2 2 2 
w max + w i ( w max + w °o) 

Thus compared to ( 38 ) only the expressions for tj_ and tg are exchanged. 
It is easily checked that t 2 > 1 > tq is exactly fulfilled for 

> w^, and that for > w M actually t-^ + tg - 2 < 0 , as it must be. 

With it the problem of defining 3, tj_, and tg for prescribed 
ratios v-j. : w M : is completely solved. For < v M equations 

( 36 ) and ( 38 ), for >v M) ( 36 ) and (40) must be applied. The case 
of = w m is technically of no interest and mathematically trivial. 

In conclusion, the quantity 5 is computed. 

y 2 +w 2 v 2 _ w . 2 

max w °° "max W 1 

2,2 2 2 
Wmax + Wmax " 
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1 + 3 \/' fc 1 t 2 - 


2 w 2 
cw max 


2 , 2 
Vmax + w± 


V 2 + W 2 W ^ + V 2 (w 2 + w 2 - 4 w.w ) + W 2 W 
1 2 mas 00 mas mas co i 1 00 1 « 


w, 


mas 


+ Vj_ 


/ 2 2\2 
l w max “ Veo ' 


tj + tg Vpa ^ 4 + y^V 5 ^ 2 + Vj g - WjVB,) + Wj^oo 2 


( w 2 +. v 2 )( v 2 - -w 2 ) 

v w mas w i ' ' mas 00 ' 


whence after an easy calculation 


(Wco ~ -Wj ) 2 

w mas “ w oo 


C^D 


This important formula of Ruden’s theory is obtainable without difficulty; 
it is valid for Wi = w^,. 


Another interesting and practically important result is the 
following: the quantities denoted by \r m) and w^ correspond 

to the normal operating condition and are now written with capital 
letters. When any desired operating conditions for a fixed diffuser 
are to be studied, the equation system (32) is used with general t Q - 
instead of ( 35 ); w* can then also be replaced by By 

limitation to w^ < w^, and taking (36) and (38) into account where 
capital letters are used at the right-hand side. 


WjiWooi W^ = to’ 


’max 


W_ c + W ± 2 W ^ 2 + W M 


2 W 


mas 


2 W. 


max 


*1 - % 


fT-bT f ° r = *3 


— 2 — for t 3 < t D % 1 


< t 2 - 1 


t o ~ ^1 
1 - ti 


for t 0 ^ 1 
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whence 


to = 


w i w max 2 + W oo 2 


W„ flT " + W ± 2 


v tj 2 
00 "xnax, 


v. 


When the ratio — gradually increases from 0, t Q increases 


w„ 


WrnR-v- 

correspondingly. The ratio -j decreases continuously so 


w 


long as t Q remains less than unity; then it rises again for t Q > 1. 
Thus the smallest increase of speed occurs at t Q = 1, that is, for 


w 


W 2 , tj 2 

"rryrr + w j 

y 2 . y 2 

"max 


In figure 15 the velocities at the end points of the nose of a specific 


Ruden inlet diffuser 


h 

w„ 


w, 


= 0.4; ■■■ ^S ax 111 = 1.2i are represented divided 


v. 


hy the flight speed, as function of — ; the higher of the two speeds 

"oo 

ls w max' 


The diffuser shows a distinctly demarcated range of favorable 


efficiency. This range is given by 


¥ i 


w i ^ w max 2 + ¥ i 2 


W “ v ~ 

"oo ”00 


w, 


max 


2 2 
+ Woo 


If 6 and are given, — can be computed by (4l). Figure 14 

W 1, W - 


in Ruden 1 s report (reference l) represents this relation graphically 
by the curves 


W, / W \ 

— i = constant in a ma ^ , 5 plane. This figure 

\ ¥ oo y T.r 2 + y.2 


is now complemented by the addition of the curves 


W 


max 


— j = constant replaces — = constant, ( — = constant 


W. 


¥ max ' + ¥ oo 


— = constant 


w i\ _ 


w. 


<»/n 


W, 


w. 


oo/O 
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replaces 


W 


max 


+ W, 


y 2 + y 2 

max co 


= constant In order to indicate that the 


upper and lower limit of the favorable efficiency range are included. 
Figure 10 makes it possible to define the limits of the favorable 


efficiency range for prescribed values of 5 and 


W, 


max 

Woo 


The second diffuser type for fixed w± : w m is mathematically 

very simple. The vertex of the pair of straight lines is placed 
toward (-1,0), while the straight lines themselves rise at 4-5°; t Q is 

put equal to unity; w* coincides with w^.- y . By (32) 


w i : V oo : w max 


1 + 


•1 + 3 


hence 


3 = 


'Wmax ~ ,Wqo 

w„ 


\f 


~ _ v mR.y ~ v j 

1 2 ~ w - W W 
max 00 i 


(te) 

(W 


■For 6 


8 = 


v nax ~ 
w„ 


Wi t-j_ + tg 


w max 


- W4 


v, 


max 


w, 


max 


To keep 8 at a minimum, quantity 


ti + "kg 


must be reduced to 


a minimum. The geometric mean of tj_ and tg, \J tjtg is defined 
"by (43); the arithmetic mean ^ (tj_ + tg) becomes a min i m um when 

tj_ and tg are made to come together as closely as possible. As 
tg >1 > tj_ i-s to be valid, the optimum is only appr oxi mately 
attainable. For Wj_ < w^, \jtjtg >1, and the optimum is represented by 
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(w, 


ti = 1, t 2 - tjtp - •j— 


~ w±) 


max 


max " w oo 


2 2 


F 


v ± 2 


For w-j_ > w^, these values must "be exchanged, for tj_ and. ± 2 . In any 
case ' 


t*i + tr 


V, 


max ( w oo + v± 2 ) - 2v max w 1 w 00 (v i + wj + 2w ± 2 w* 2 

max w oo' 


2 

2w ± (w m 


and. 

6 = v max^°° ~ v l) 2 
2^1 - V J 


(44) 


This value is now indicated, by 5* and then compared with the values 
given by (h-l) : 


5 _ 5 w max^ w max + v °° ^ 
* 


2wi 


w 


00 


The quotient 


is in any case greater than unity. The Becond 


diffuser type is, therefore, inferior to the first. For very Fmw.l 1 
and very great w^ : w^ &* must become substantially greater than 5. 
Because for w. < w : 

1 ^ 00 


w, 


( v may w oo) 


max 11 "max 


2W ± W*, 


oo/ Wo^Woo 'W’oo 

> 2w 1 Woo “ W ± 


(45) 


and for Wj_ > w 0 


^Trw.x^mex + v ±( v ± + v oo) 


2Wj W- 
1 00 


2^1 w =o 



( 46 ) 
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so that — constitutes a measure for the minimum increase which S 

w_ 

00 

experiences with respect to 5. In (45) and. (46), w^y was replaced 

"by the greater of the speeds Wj_ and w M . Consequently, the estimates 

are much "better as the speed increases are smaller. But small 
speed Increases correspond to small parameter values 0. 


4.- Theory of Optimum Diffusers for Any Variable Ratio Wj_ : w^ 

In thi s instance three types of diffusers are involved. For the 
first two, C has a straight central piece that lies on gg or g-j_; 

for the third type the central piece of C is curved and touches gq 
a.nrl gg . The straight line g runs parallel to the t-axis, while gg 

passes through the origin of the coordinates . 

The first type was theoretically discussed; it "belongs to the Ruden 
group. Now, however, the reference point t D must "be chosen general, 

R.nrl at the same time the relation tg =0, that is 

2V2 - h - * 0 


must "be observed. By (34) 


2 4(tj_ + tg “ 1 ” tj_tg) 

(t 2 - t ± ) 2 

the insertion of 



according to (47) gives 




3 ^ 


BACA. 1M 1267 


or differ ently expressed 


P 2 = 


t, t 
1 2 


Accordingly- 8 is computed as follows: 


8 = 


1 + 


_j ri + t 2 

2 




Since w^ 


w the velocity ratio follows as 


Wi 




for t D 


for t Q 


<1 
£ 1 


Being chiefly interested in the case of very high values of w 1 : 
t Q is accordingly put as t G — >00, so that 


v max • 

It is then readily apparent that 
Therefore 



tj_ must be greater than i. 



However, since speed increases of this order of magnitude are 
inadmissable in practice, the first type is unsuitable for the present 
purpose . 
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hence 


For the second type, f3 = 1 and tj_ + tg = 2. Again = w*, 


- 1. 


V • TjT • V 

i 00 ‘ max 




1 + \ / 1-, t, 


1 , 

2 : < t, 


t2 “ 1 


for t Q < 1 


'12 


o t l 


for t 0 s’ 1 


For very high t Q values 


1 + \7 t x t g 1 + \/T - a g 

w max : w i ~ 1 _ + “ a 


with t^_ = 1 - a, tg = 1 + a. By proper selection of a the last 
ratio can he made to approach unity as closely as desired, so that 
in this respect the conditions are much more propitious than in the 
former case. On the other hand, the conditions for medium speed 
ratios are still too unfavorable. For from 


t o ~ 


1 + nAi^ V 1 1 +\/l - a 2 V 1 


w m 

00 




v i 

when assuming — < 1 and hence t Q < 1, follows 


w max - 2 


t 2 ~ t c 
*2 - 1 


= 2 


1 + a - t_ 




= 2w_ + —I w„ 


1 +n/i~ 


w ± > 2w a 


a result useless in practice. 

The third type is characterized by the fact that C touches 
g^ and gg 1 . The points of contact have the abscissas -t^_ and -tg. 

The equation of the central arc of C is by ( 29 ): 
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s 2 = (t + l) 2 - 0 2 (t + ti)(t + trj) 

which, after differentiation and putting t = -t]_, gives 

2s = 2(1 - t-L) - P 2 (t 2 - tp) 

t=-tp 

By assumption this expression must disappear, whence 

P 2 = 2 ( 48 ) 
t2 " tj_ 

This value, introduced in the expression of s 2 and differentiated 
for t = -t 2 , gives 

2a ff = 2(1 - tg) + 2 l - %■ (t 2 - t ± ) = 2(2 - tj. - t 2 ) 
t=-t 2 t 2 ' *1 

or, since s(-t 2 ) = t 2 - 1 j 

2 - tj_ - t 2 - 

t2 ' 1 

In order that C touch g 2 this value must he equal to 
t2 - 1 

that is, 

t 2 

(2 - t^ - t 2 ) t2 + (t 2 - l) 2 =0 
or 



t x t2 - 1 


(^ 9 ) 
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By (48) and (49) 


P 2 


2t l 
1 + t! 


or 



§L_ 

2 - p 2 


Quantity 6 is expressed by P in the following manne r : 


5 - —P / P 2 

■ 2(1 + P)U . p 2 


+ g..--.,..P 2 

P 2 



2(1 + P)(l - g) 2 
P(2 - P 2 ) 


The connection between P and 5 is numerically represented in 
the table as follows: 


JL. 

1.0 

0.9 

0.8 

0.7 

0.6 

0.5 

o.4 

5 

0 

0.03 

0.13 

0.29 

0.50 

0.86 

1.39 


To keep 5 within tolerable limits , P certainly must not fall 
below 0.5. Formulas ( 32 ) then give (w^ = Wmp Y ) : 


t Q 

v i : v oo •' v max = “ 


2 - P 2 - P 2 t c 
2(1 - P 2 ) 


for t Q ^ 1 


1 + p * 1 + p 


< 


(2 - 02) t Q - P 2 
2(1 - P 2 ) 


for t Q = 1 


that is 
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and. 


Va - - g(i . p) w i for V 1 S v « 


_ 2 ~ 3 s 


£. 


> ( 52 ) 


w max = 2(1- p) w i ‘ 2 ~(l - P)' W “ for Vi ^ W °° 


for P = i, for example, 


v max = < 


1.75w M ~ 0-25vi(w i -g w a ) 
l-T5Wi - 0-25w 0O (w 1 > w^) 


The conditions are seen to he much more favorable than in the other 
two cases. Nevertheless it is desirable to reduce the speed increases 
still farther without increasing 5. This is accomplished, geometrically 
speaking, by placing the central arc of C a little lower. Analytically, 
this operation signifies a decrease of p 

and tg. The corresponding wall thickness 

decreases with decreasing p, and, ( 51 ) being valid, the velocity w* 
follows as 


for fixed values of t n 


S = 


jHr P\ 


'^i + ^2 


v* = 


= 1 + 3 

1 - t- 


- (w - tiWi) for w, < w 


w. 


- ^ ( w i " t l v «) for w i = 


> 




(53) 


A decrease of P therefore also acts favorably on w*. For extreme 
values of w^ : w^ Hy = w* is, of course, no longer valid, and the 


range of validity of the last equation diminished with decreasing p. 
The greatest increase of speed grows, in any case, when p decreases. 
The extent to which P is to be reduced is a matter that mudt be 
decided in each case individually. 
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III - EXAMPLE TO II, 3 


It concerns Euden's diffuser II (reference 2) with constant 
increase of speed =. 1.2 along the nose contour in normal 

operating condition = 0.4. By ( 36 ) and (33), la this case 3 = ^ = 0 . 18 , 


ti = = 0.807 and 

nose is represented by 


t 2 = 


122 ' 


= 24.4, so that the upper contour of the 


X = T - log T 


y = 0.18 \/(t " 0.807)(24.4 - t) + 4.55 arc tan V T - 0.8 01 

\f 24.4 - t 


- 1.6 arc tan 


11 Vt - 0.807 

2\]2k.k - t 


h = 1 . 8ir = 5 . 655 = one -half inside width . 


Figures 11 to l4 give a comparison of the computed and measured 

pressure distributions for four different — ratios. The ordinate 

is plotted against the abscissa x, (as in reference 2 ) that is, 
the ratio of static pressure p - p M to kinetic energy at infinity, 

q = £«• 2 . By Bernoulli’s equation 

VL _ >** 


^st 

q. 


p + £ w 2 = ? +■£■ v 2 

= 2 2 05 


that is, 
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The quotient — is defined "by the graphical method represented in I, 5- 

^co 

The agreement between theory and test Is good. The larger discrepancies 
at the diffuser exit are explainable by the effect of the blunt end, 
while defects in workmanship are certainly noticeable at the mouth of 
the diffuser. The only appreciable point of difference is that the 
marked low-pressure peaks at' the mouth are blunted by friction effects 

as postulated by theory for very high and very low — . The exact value 

Woo 


of the maximum low pressure was, of course, not measured, but the 
measurements give a fairly good idea of the pressure distribution along 
the mouth of the diffuser.. According to theory the maximum low pressure, 
that is, the maximum velocity of the flow, lies always at one of the 
end points of the nose contour. The velocities at these points, 
each divided by the flight speed, are represented in figure 15 


for variable — ; is the greater of "the two. The favorable 

W oo 

effective range, filling the interval from 0.4 to 0.66, is plainly 

■Vf—o -y- 

visible. The dashed line represents the ratio for w-» > w„,. 

w i 

Continued toward the right it reaches the height of 9 . 3 . 


IT - EXAMPLE TO II, 4 


w ± 

Three different types of diffusers for unlimitedly variable — 

CO 

were discussed, but the first two were rejected because of excessively 
high increase of speed even with great wall thickness. The third type 
depends only on one parameter 3. The values t^ and tg for given 3 
are defined by (48) and (49). But, as stated previous'ly after t^ 
and tg are defined, 3 can be subsequently made variable again, 

so as to provide more favorable conditions for the operating range of 


principal interest by foregoing extreme values of — . 

w„ 

00 

start is made from 3 = 0.5, for which t^_ = ^ and 1 2 

ti +..tp has the value 7.14. This value is replaced by 
matically convenient value 


Suppose the 

- 7- The sum 
the rnathe- 


(54) 


tj_ + tp = 8 
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which is equivalent to a slight decrease of the initial value for P; 
this is now 0.47. With 


tq = h - \/l5~ = 0.127 


and 


tg = h + \/l5 = 7.873 


to "be computed from ( 5*0 and (49) the equations ( 17 ) for computing 
the diffuser contours are set up: 


1 + 8 arc tan — . 1 . P -Ur Z, - 2 arc tan 7.873( T ~ 0.12 Xl 

V&r - r2 - 1 V 8r - te - 1 


In accord with previous studies, the values chosen for P are 
the initial value O.hj. Figure l 6 shows the contours for P = 

0.2, 0.3, and. 0.4. Scale variations ensure that all diffusers 
equal absolute wall thickness d. The values of the relative 

thickness 5 = % are also indicated, 
h 

Figure 17 shows the maximum speed distribution for the four 

diffusers plotted against — , which is, for w., < w'oo the ratio 

w 1 

CO 

, for w-t > w„. the ratio . A decrease of P in the 

w ra J i a w ± 

central range is favorable, at the ends, unfavorable. The others 
refer to the case P = 0.4. Figures 18 to 20 represent the relative 
pressure distributions for different operating conditions in comparison 
with the pressure -distribution curves of figures 11 to 14. 


below 

0 . 1 , 

have 


x = t - log r 


y = P 


\J8t - t 2 - 


h = (1 + P)jt 
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V - APPENDIX 


1.- Variation of the Function ~ of t 


Writing x, Xp, xp, and. y instead of t, tp, tg, and ^ 

so as to conform with the conventional notations of analytical geometry, 
the equation of the central arc of C reads 


y 2 = (x + l) 2 - 0 2 (x + Xj_)(x + Xg) 


x 2 (l - 3 2 ) - y 2 + x 2 - P 2 (x-j_ + Xg) + 1 - p 2 x-|_Xg = 0 


If a curve of the second order 


O p 

a^jx + 2a 12 xy + aggy^ + 2a 1 jX + 2a 2 ^y + a^ = 0 


is to he analyzed, the determinants 


an 

a 12 

a 13 


a 21 

a 22 

a 23 

(ait = a ki ) 

a 3l 

a 32 

a 33 




a ll 

a 12 

A 33 

= 

a 21 

a 22 


must he home in mind. The decomposition of the curve in a pair 
of straight lines is indicated hy A = 0, while for A ^ 0 
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^3 


■^33 > ^ gives an ellipse 


=0 a parabola 


A = 


1 - P c 


A^^ < 0 a Hyperbola. In the present case 

Zp+Xgl 


1 - P c 


-1 


2 

0 


! - P 2 0 


1 - P 2 XpX g 


5 2 X 1 + V 

2 ) 


= ^1 - P" 

- (1 - P2)(l - fS^pXg) 
= p 2 (l + xpXg * X 1 “ x 2 ) 

+ ■ x i' £ 


*2 - ^ll 


P 2 (P ? - Po 2 ) 


with 


_ 2^(1 - x 1 )(x g - 1) 


(x g - x 


o7 


Decomposition therefore occurs for P = 0 and. P = P Q , where, 

since 1 - P 2 = — X j- - + ^ — 

(xg - x 

and attains unity only for xp + Xg = 2; A33 ha.s the value P 2 - 1, 

hence ^ hyperbolas when P < 1, ellipses when P > 1. The case P = 1 
results in a parabola, or, for Xp + Xg = 2 in a pair of parabolas, 

of which only the one straight line above the abscissa axis is in 
evidence. 


2 ) 2 

D 5 " 


is certainly not greater than unity 


2 - General Determination of Wj 


Elsewhere w max had been defined geometrically. For the case 
that Wmax is assumed at an inside point of the nose contour, 
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a tangent is placed from a given point- of the t-axis on the conical 
section to which the central arc of C relates. The corresponding 
calculation is as follows: The equation of the conical section is 


(t + l)^ - 3 2 (t + tj_)(t + l^) - s2 = 0 


or 


(1 - p 2 ) t 2 + 


[- 


P 2 ( t-i + to ) 


t + 1 - 0 2 tq_t 2 “ s2 as 0 


while the considered point has the coordinates t = -t Q , s = 0. 

For computing the tangents from point (x Q , y Q ) on the conical section 


au*2 + 2a l2 xy + a^gy 2 + 2a^ 3 x + 2a g3 y + a^ = 0 
the "polar" of (x Q , y 0 ) with respect to the conical section 


auxx 0 + a 12 (xy o + + a 22 yy Q + a 13 (x + Xq) + ag 3 (y + y D ) + a 33 = 0 
are determined. 

This straight line dissects the two constant points on the 
conical section. In the present case the equation of .the polar 
is 


-(1 - P 2 ) tt Q + 


1 


- P 2 


ti + tr> 
2 


(t - t Q ) + (1 - P 2 t-jtg) = 0 


(l - 0 2 ^ ■ * tg ) t 0 - (1 - P 2 tqtg) 
(1 - P 2 ) t Q - ^1 - P 2 ^ g ^ j 


( 55 ) 
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The polar is normal to the t-axis . Denoting the value given 
by (55) with -t 0 ', gives 


v max 



(l + P) 


that is 


v max ) / b «.2 (tp - tp') 2 

sW) 


A simple calculation gives: 


'v \2 

"maxi 


Mi + P) 2 


Ho 


(t2 - tj 2 P 2 (P 0 2 


[■ 


-5-' (1 - to) 2 - 3 2 (t x - t 0 )(t 2 - t Q ) 

r) 


(56) 


P Q is the value given by (20); 3 is smaller than 3 0 . The formula 

holds for the case that the value tg ‘ ranges "between t^ and t 2 , 
hence for 


t! = 



t Q " ) 

i 1 - e 2 V2) 

(1 - p2) t Q - 

(*- 

p? M + t 2 N ' 
2 J 



= . t 2 


Determine next the two t Q values for which t Q 


t Q 1 = tg. These values are given by 


£ 


t 3 = 1 


7T (tg - t x )(l - t x ) 


£ 


1 ” t^_ - 2 (tg ” tj_) 


and 


£ 


% = 1 - 


(t2 “ t x )(t2 - 1) 


P 2 , 

t2 - 1 - — (tg - ti) 


= tn or 


(Figure 2l) 
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For 


P 2 < 


2(1 - t x ) 

t 2 - t;L 


ranges •'between 1 and °° , for 


2 g (*2 - X > 

P < t 2 (t2 - tq) 


t^ lies between 0 and 1 . If 3 satisfies these two inequalities, 
formula (5 6) holds for all t Q that meet the conditions 0 ^ t Q ^ t^ 
or t Q ^ t3, while w^.y = w^ for t^. < to g t^. If only one of 

the inequalities for 0 is satisfied, there is only one validity 
interval for (56) , and if none of .the inequalities applies, Wm^ = w # . 
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PART II 

I - AUXILIARY MATHEMATICAL EXAMINATION 


Consider the conformal mapping effected by w = z + — 

z 

in such a way that z and w are made to vary in the same plane, 
z being restricted to values which, geometrically speaking, lie 
outside of the unit circle and above the real axis. The case then 
presents a conformal mapping of a semicircularly notched half plane 
on the entire half plane. 


It is seen at once that 

(a) each inside point is shifted nearer to the real axis 
by the mapping 

(b) the boundary points lying on the real axiB all travel 

outward, that is, for z > 1 toward the right, for z < -1 
toward the left, • 

('c ) every boundary interval lying on the real axis is reduced 
in its length. 


These three facts hold, as will be shown, not only for the special 

mapping according to z + i, but very generally for the conformal 

z 

mapping on a half plane provided with any desired notch on the full 
half plane. More exactly the following is valid: (fig. 22) 


Let B represent a simple connected space in the upper half of the 
z -plane bounded by the two^ semi straight lines z < a -and z ^ b 
(a and b real, a < b) and a curved arc free of double points 
connecting the points a and b and with exception of the end 
points running entirely within the upper half plane; B is mapped 
by the function w = f(z) on the whole upper half plane and the 
mapping function existing according to Riemann's mapping principle is 
so standardized that the infinitely remote point transforms in the 
finitely remote point and that in the development 


Cn C q 

v = Cz *+ c Q + — + — — + . . . (C > 0) applicable to the vicinity of 
Z z 2 

infinite distance, the coefficients C and c 0 attain the values .1 
and 0 (standardization at rest in infinity, w = z + ((0)). 
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In. this instance 

(a) J(w) < J(.z) for all inside points 
(4) w > z for z ;> b, w < z for z < a 

(c) w 2 " W 1 < z 2 “ Z 1 f° r z o > Z 1 ^ 4 and for z, < zo < a 
(W! = f(z^)) 1 x 

Proof: 

For the first it is noted that the function w(z) can he 
analytically continued by the reflection principle beyond the 
real axis. Putting z = x + iy. 


w - z = w' - u' + iv ' 


and considering v ' as function of x and y, v ' is a regular 
potential function in the entire space B, which on approaching the 
infinitely* remote point disappears . Along the boundary piece on 
the real axis v* = 0 , along the rest of the boundary 
v' < 0. Because of the regular behavior at infinity the validity 
of v ' < 0 on the inside can be deduced from the ine quality v ' < 0 
at the boundary. Since v' = J(w) - j(z) the claim (a) is therefore 
proved. 

Let (x 0 , 0) represent a boundary point of B on the real axis. 
As v^Xqj 0) = 0 and v’(x 0j y) < 0 for sufficiently 

small positive y values (x Q , 0) is certainly < 0. 

The equal sign must be excluded, for if (xq, 0) were equal 

to 0, then as well as would vanish in point (x 0 , 0), that 

is, the derivative would be 0 at the particular place and the. 

development of w’ would assume the form 


w' = ccq + a K (z - x 0 ) K + . . . 

with 


K >2, a K ± 0 (a Q real) 
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Bat then to a small two-dimensional surrounding of x q there would . 
correspond a multi -lohed surrounding of point ccq through the 
mapping according to w'(z), and v 1 would have to assume positive 
values in the upper half plane, which cannot "be, according to (a). 

By the Cauchy -Riemann differential equations Tp— is <0 for all 

boundary .points of B lying on the real axis , that is, u 1 
decreases when traveling from left to right. As to the 
disappearance of u' at infinity it is seen that u' must be positive 
for z ^ b, whereas u' assumes negative values for z < a. But 

along the real axis u* = w - z, whence the claim (b) . 

Lastly, for z 2 > Zj_ ^ "b and for z i < z 2 ^ a: 

(w 2 - Wq) - (z 2 - Zq) = (w 2 " Z q ) - (wq - Zq) = V-'g ~ u'q 
This monotonic behavior of u' proves the correctness of c. 

II - MATHEMATICAL TREATMENT OF GENERAL INLET DIFFUSERS 


The following investigation contains four complex variables 
all of which are in analytical relationship to each other. The first 
variable, z , is the complex coordinate in the plane of the diffuser, 
the second, W, is given by the value of the stream function, the 
third, w, by the value of the complex velocity, and the fourth, t, 
is a pure mathematically explained auxiliary quantity. The range of 
variation of z is represented in figure 23(a). The diffuser is 
assumed to be of infinite length; it is to be symmetrical and bounded 
by two convex curves with constant direction of tangents, whose pieces 
extending to infinity are rectilinear from a certain point on. The 
convexity required for the outside is natural; but at the inside wall 
the limitation to convex formB imposes a relingulshmfcnt not shown 
before ha nd by physical-technical considerations. On the .contrary, 
the constructed engine cowlings used in practice for wq < rather 

exhibit cross-sectional enlargement downstream from the diffuser opening, 
that is, no convex contours (fig. 23(b)). Nevertheless the limitation 
is restricted to convex forms for the reason that the mathematical 

analysis affords a simple connection between quantity S = ^ (figs. 23(a) 

a.nri 23(b)) -and the velocity distribution along the contour; but the 
quantity 8 is decisive only for convex forms in problems of engine 
installation, while on forms with cross-sectional enlargement the 

quantity 8* = — is decisive the dependence of which on the velocity 
h* 

distribution is far more complicated. 
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, ¥i 

Only on© flow which corresponds to a specific ratio — =■ is 

investigated. It is represented by an analytical function ¥ = F(z). 
Putting ¥ = IT + iV, the streamlines are given by V = constant; 

w = is the complex velocity whose conjugate complex value 


represents the actual velocity in magnitude and direction. The 
function F(z) can be replaced by cW + c ' with positive c and 
any complex c f without altering the flow pattern, because only 
a multiplication of the complex velocity by a constant positive 
factor is effected, which is equivalent to a variation in mass unit. 


Ely ¥ = F(z) the z space is mapped on a symmetrical slot region 
which is, the slots run parallel to the real axis 

rightward to infinity. By taking advantage of the previously cited 
freedom in the determination of ¥, it can be assumed that the two 
slots are given by V = ±it, U > 0 (fig. 2*0 . 

The quantity w = varies in a symmetrical Bpace lying 

dz 

entirely within finite limits. Putting w = pe^ T , the actual velocity is 
indicated by pe~i T . The boundary values of t are given directly by 
the diffuser contours, -t = arg dz, where. It is true, the orientation 
of the line element still remains questionable, that is, t is defined 
up to multiples of jt. In any case it follows from the assumed convexity 
of the contours that t varies monotonically, if the contour Is 
followed from the stagnation point to infinity in one or the other 
direction. The further study is restricted to the upper half of 
the z space. The image of this half space in the w -plane is bounded 
by a length lying on. the real axis, which corresponds to the 
line of symmetry of the diffuser and the two straight walls, 
and by a curve arising from the curved part of the diffuser contour. 

The curve proceeds entirely in the upper w-plane when the stagnation 
point lies on the outside diffuser wall; It runs In the lower half plane, 
when the stagnation point lies on the Inside wall; It splits into two 
arcs, each belonging to a half plane, when the stagnation point Is 
a point of the curved contour. In the latter case the two branches 
of the curves meet in the zero point of the w-plane and have a common 
tangent for w = 0. In any event the curved part of the boundary of 
the w-space is free from double points because of the fttonotonic 
variation of t ; from this it can be concluded without difficulty 
that the w-space must be a smooth, that Is, single lobed, space. The 
figure 25 represents potential space forms. The zero point Is always a 
boundary point. The w^ and w» indicate the velocity Inside the 
diffuser and the flight speed. 
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When the w-space with the points w* and Vo, is given, the 

diffuser design is defined. For the w-3pace is mapped on the W-space 
of which only, the upper half is taken into account; w = 0 changes 
to W = iti, while to wp and w^ correspond the two infinitely 

distant "boundary points of the w-space. The invariably existing 
mapping function W(w) enables z to be computed as function of w. 

Since 


w 


dW 

dz 


z 



P 1 dW . 

= / - r— dw 
/ w dw 


According to "it z appears to be determined up to an additive constant. 

The other possible similarity mappings in. the z -plane follow when it 
is remembered that the quantity W is .arbitrarily normalized by the 
foregoing, and that truly directional similarity mapping with 0 as 
central point in the w-plane is admissible. 

Since pressure and velocity distribution are the chief factors 
governing the quality of the diffuser, the usual process is to 
proceed from the velocity patterns and to determine the contour 
according to it. 

In the present report the hodographic method is resorted to. 

But, while in reference (l) the amount of the w-spaqe was automatically 
restricted by the assumption that the curved part of the contour is 
a half circle, no special assumption is made here. 

It is convenient to replace the W-space by a half plane, putting 
W = t - log t + jt i - 1. It is readily verified that by this relation 
between t and W the W-space is, by appropriate fixing of the logarithm, 
mapped on the upper half of the t-plane. The two infinitely distant 
points of the W-space correspond to the points t = 0 and t = ^ 
while W = iti changes to t = 1 (fig. 2 6 ) . 
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The w-space must accordingly he napped on the upper half plane 
hy means of a function t(w), the points 0, v±, w M must change 
into 1, 0, co, (fundamental mapping of w-space with boundary normalizing). 
The corresponding diffuser form is then defined hy 


z = flld w = r~( 1 - A — dw 

J w dw J w \ tydw 


III - A GENERAL THICKNESS FORMULA 


The main object of the present note is to prove the fact that 
for fixed "values of w^, w^ and w MT the corresponding wall thickness 
of the diffuser cannot become less than the value computed by Ruden 
for the case of constant increase of velocity along the nose 


(Vqo - ~Wj ) 2 
2 2 
v max " w °° 

only Ruden' s diffusers reach the minimum value. 

To be able tOs compute the quantity & for any desired w-space 
only the behavior of the function z(w) near v M needs to be known. 

The quantity t becomes infinite for w M in first order and has a 
development of the form 



+ reg. 


with f positive for "»! < w^ and negative for w^ > w^. It is 

to be noted that t(w) over an interval of the real apcis 
containing the point w^ on the inside certainly can be continued by 
the reflection principle. 
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Accordingly 


t 


•w - v„ 


(1 


(( 0 ))) 


1 

t 

dt 

dw 



(1 + 


(( 0 ))) 



+ reg. 


l\ M 
t ) dw 




+ reg. 


where reg. indicates a function regular for w^ ( ( 0) ) one regular for 
w ra and disappearing in w^ itself. Adding the development of — 


1 


w 


w 


(w 




W J + • • • 


gives then 


z 



W ) 
00 



+ 


1 

W - V 

00 


+ reg 


dw 


1 


w w - w 

CO CO 


log (w - wj 

w 2 w ” j 


+ reg. 


From this formula the amount of d + h, that is, the variation of the 
imaginary component of z on passing through w^ can he taken at once 




d + h = it 
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Quantity h is correspondingly given "by the variation of the imaginary 
component of z on passing through wj. In this instance 


t = c(w - Wj_) + . . . c f 0 


111 


t C W - 


+ reg. 


! _ 1\ dt 


tj dw v - Wj_ 


+ reg. 


I = i- + ((0)) 


V w. 


z = 


± + ((0))l 

W 4 


w - w^ 


+ reg. 


dw = - — log (v - w. ) + reg. 

v., 1 


that is, h = — : hence 
v i 


d + h 


8 + 1 = (f + wj 
w 

00 


and 


8 



+ WjW 
-L 00 



While in the. formulas for h and d + h the arbitrary standardization 
of W plays a part, the formula for 8 is of completely general 
validity. The quantity f, that is, the residuum of t in point w ro 
has the dimension of a velocity. At fixed values of w^ and w^ °° 

8 1b solely dependent on f. 

The thickness formula is now used to compute the wall thickness 
of a Ruden type of inlet diffuser. The w-space is a half circle; the 
center is the zero point, the radius is and for w-j_ <w ro the 

half circle lies in the lower, for Wj, > w in the upper half plane. 

The mapping function t(w) is a rational function of the second degree, 

v 2 

with zero points at Wf and while assuming the value » for w M 

T.r . 
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v n 


and. 


ma,x 

w_ 


Accordingly 


t = constant 


(w 


y iH g : — - 1 


(v - wj 




and the constant factor gives the value 1, since t(0) = 1. The 


residuum for w„ is 

CO 


f = 


(v - v ± )lw - 


V, 


im.Y 


00 Wj_ 


Vmax 

W co “ W„ 


= ^00 “ W i) 


("max 2 - vvi) 


("max 2 - w 2 ) w ± 


hence 


5 = i. 

w 


2 _ 


, x "max - ■w e ^ 1 

( w oo - v i) 7T + v i " 


V 2 - w 2 
w max 00 


_ (voo - w ± )‘ 


"max' 


2 - v 


IT - EROOF OF THE MINIMUM FORMULA 


A piece is tut oup of the w-space B of figure 25(a) (w^ <v m 

stagnation point at inside -wall), as indicated in figure 27(a). What 
is the variation experienced ty 5? According to the cited thickness 
formula it is sufficient to analyze the variation of f. It is f >0. 


The reduced space is denoted ty B ' . The function t(w) that 
maps the original space on the half plane is regarded as known. The 
mapping function t'(v) for B' is then obtained in the following 
manner: first map B' ty means of t(w) ; the result is a half plane 

with a notch. This is then mapped ty means of a function t"(t) on 
a complete half plane, the mapping to te standardized at infinity 
according to the formula t" = t + ((0) ); t 1 is then a whole linear 
formation of t", t' = It" + p. Hence, since t' as function of t 


at infinity has a development t' 
f r = Xf holds for the residues f 


= .It + p + + —■ . . . , the relation 

and f '. 


t t 2 
Quantity 


X is defined ty 
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the condition that points 0 and Wj_ shall change to t ' = 1 ' and 

t’ = 0. Hence t(0) = 1 and t(w^) = 0. At transition from t to t 
the interval (0,l) is reduced in any event, since the notch lies outside 
the interval. Accordingly X must definitely he greater than unity, 
hence 


f ' > f 


The wall thickness therefore increases as the w-space is reduced. 

When the w-space .increases, 6 decreases accordingly. So, if a specific 
w max ts Siven, which may he attained hut not exceeded, the w-space 
may not extend heyohd the circular space | w j ^ w mT . The wall 

thickness is a minimum when the space available is completely utilized, 
that is, when the amount of the velocity along the entire nose is made 
equal to w^. This proves the minimum theorem for w 1 < v w when 

assuming the stagnation point located at the inside wall. 

Next, the stagnation point is placed on the curved portion of 
the contour, while w^ < w^ (fig. 27(h)). In this instance any increase 

of the w-space in the lower half plane is again associated with a 
decrease in 6. This decrease in 8 can he continued further hy 
reducing the part of the w-space situated in the upper half plane. 

This is readily verified hy the mapping t'(w) of the reduced space 
as above, passing through t(w) to a half plane with a notch, then 
completing the transition to t" = t + ((0)) and to t' = Xt + \i. 

The difference now consists in the notch in the t-plane lying 
between 0 and 1, thus increasing the interval (0, l) to t" on ■ 
transition. Therefore X < 1 and f 1 < f , that is 8 becomes smaller 
as claimed. The entire piece of the w-space lying in the upper h a l f 
plane is omitted, wherewith the stagnation point shifts to the inside wall. 

The optimum within the group of convex diffusers for specific w^-y 
is therefore actually reached hy Euden's inlet diffuser. 

If > w M the proof is entirely analogous. The residuum f 

is then negative. A notch in the w-space of the upper half plane appears 
then in the t-plane between 0 and 1 and results in X < 1. 

Accordingly f = Xf > f , and 8 increases. It Is seen further 
that a piece of the w-space lying in the lower half plane must he 
reduced, in order that 8 decrease . The optimum w-space for a given 
’ w max lB therefore the half circle |w| J(w) >0. 
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V. CONCLUDING REMARKS 


The subsequent discussion is limited to the case w^ < w M . 

As already stated, the diffusers used in practice exhibit cross- 
sectional enlargements downstream from the diffuser orifice (fig. 23(b)). 
The corresponding hodograph record is given in figure 28, it being 
assumed that the stagnation point lies on the In side wall. The 

quantity 8 = — can be reduced below Ruden's value, if the holographic 
h 

picture is so chosen that, aside from the straight piece, the lower 
half of the circle | w | = w and between 0 and wj_ an arc 

running in the lower half plane takes part on the demarcation. 

The nearer the last arc is pulled to the remaining part of the 
rim, the smaller 8 ■ becomes. The value 8 = 0 is attainable. However, 
two facts should be borne in mind. 

First, the quantity governing the practicability of the diffuser 
d* d 

Is 8* = — , not 8 = — , so this quantity would have to be analyzed In 

relation to the w-space. The resulting relation might not be as simple 
as for 8. In any case it is doubtful whether 8* varies in general, 
as 8, and to what value 8 my be reduced. 

Second, the study of inlet diffusers of constant Internal 
cross section can be restricted to Wj_ : w w : w mT and 8, but not 

on "the forms considered here. Ruden's diffusers indicate considerable 
regulating capacity and operate nearly without loss to the extent ' 
that the conversion of velocity to pressure Is essentially effected 
upstream from the diffuser mouth. To forestall 'excessive variations the 
section enlargement must be done very carefully, that is, the holographic 
pattern must not depart very much from Ruden's. In any case, constant 
increase of speed along the curved part of the outside contour must be 
specified for the design operating condition, and the smallness of the 
notch ft (fig. 28) vouches for the _ outside contour itself not 
deviating perceptably from that by Ruden. To this extent the minimum 
formula proved in this report Is of practical significance in spite of 
the limitation to constant cross sections. Forms with cross-sectional 
enlargments are to be discussed in a subsequent report. 

Translated by J. Yanler 
National Advisory Committee 
for Aeronautics 
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Figure 1.- Schematics! representation of the diffuser contours and 

corresponding flow. 



Figure 2.- Schematical representation of a flow in the parallel strip, 
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Figure 3.- By dz/d? the parallel strip is mapped to a region 
symmetrical with respect to the real axis 'bounded by an 
ellipse, a finite and an infinite piece of the real axis. 



Figure 4.- By t = e" 5 the parallel strip is mapped to a plane 
that is cut along the negative half of the real axis. The two 
boundary pieces corresponding to the interval (-tg, -tj) will 
later yield the curved parts of the diffuser contour. 



Figure 5.- | dz/d£ j plotted against t(t < 0) for various values of 

P (ti_ and tg fixed) . 
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Figure 6.- The zone of variation of. the complex velocity w 

to P. Ruden. 


according 




Figure 7.- The conformal mapping of the t-space on the w -space. 
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Figure 9a, b.- If w max is to stay beneath a certain limit, the curve C 
must not approach the t-axis too closely. For fixed values of Wj_ 
and w^a) , as well as for variable values of w^ and w co (b), there 

always results a space bounded by three straight lines from which C 
must not emerge. 
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Figure 10.- The favorable effective range of a Ruden inlet diffuser is 
given by ( . The figure represents the two 

\ Woo / w 00 \y^03 

boundaries as functions of w max /w' cc and 5. The curves 
( w i/ w ») 0 = const, show in the neighborhood of w m ax/ w „ = 1 
a variation similar to the curves (w^/w^)^ = const. They are not 


fully drawn for reasons of Clarity. 
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Figure 11 


Figures 11-14.- Pressure -distribution curves for a Ruden inlet 
diffuser. Normal operating condition: w^/w^ =0.4; 

w ma x/ w ® = T ^- e s0 ^ curves represent the theoretical 

pressure distribution along the outer wall’, the dashed curves 
the theoretical pressure distribution along the inner wall. 

The circles correspond to the measured values. 
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Figure 12. 
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Figure 13 
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Figure 14. 
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Figure 17.- Maximum velocities for the diffusers shown in figure 16. 
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Figures 18-20.- Pressure-distribution curves for the inlet diffuser 

p = 0.4 of figure 16. 







Figure 21.- For calculation of w 












